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Abstract 

Election  interactions  with  interface  phonon  moaes  ana  strictly-confined 
bulk  phonon  modes  are  considered  for  the  first  time  to  calculate  the  resonant 
magnetopolaronic  is-2p+  transition  energy  of  a  hydrogenic  impurity  in  the 
quantum  well  of  a  double  heterostructure.  An  interaction  gap  is  predicted 
between  the  bulk  LO  and  TO  frequencies,  in  contrast  to  the  Frohlich-type 
interaction.  These  results  are  in  good  agreement  with  recent  experimental 
data . 


1990  PACS  Numbers:  71.38,+i,  73.20.Dx,  71.55.Eq,  63.20.Kr 
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The  interface  and  confined  bulk  phonon  modes  in  heterostructures  and 

•  •) 
sperlattices  have  been  discussed  theoretically'  and  observed  experimentally*' 

in  the  past.  Onlv  recently,  optical  phonon  modes  supported  bv  a  semiconductor 

double  heterostructure  . DBS)  have  been  solved  independently  in  the  long 

3  4 

wavelength  limit  by  two  groups.  ’  It  is  found  that  there  exist  two  types  of 
phonon  modes,  confined  bulk  modes  and  interface  modes.  The  confined  modes  can 
be  either  longitudinal  optical  (LO)  or  transverse  optical  (TO)  with 
frequencies  u>.  and  identical  to  those  of  the  bulk  excitations,  while  the 

L*  K 

interface  (.IN)  modes  can  have  frequencies  between  u>,  and  eigenvectors  and 

dispersion  relations  for  ail  these  mooes  can  be  found  in  Refs,  i  and  4.  The 
Hamiltonian  operator  describing  the  electron-optical- phonon  interaction  has 

also  been  derived  to  study  the  polaron  mobility  and  magnetophonon  resonance 

4  5 

spectra.  and  the  polaronic  states  in  a  DHS .  The  importance  of  interface 

modes  are  clearly  demonstrated  in  these  calculations. 

The  resonant  magne topo laron  effect  on  the  ls-2p+  transition  energy  is 

found  to  be  well  below  in  far- infrared  photoconductivity  measurements ^ 

carried  out  for  a  donor  impurity  doped  at  the  center  of  GaAs  quantum  wells  in 

a  GaAs/AlGaAs  multiple-quantum-well  (MQW)  structure,  in  contrast  to  the  bulk 

case  which  has  been  quantitatively  accounted  for.'7  As  the  electron  does  not 

couple  to  TO  phonons,  the  data  have  been  regarded  as  a  significant  deviation 

from  reasonable  expectations  based  on  the  Frohlich  model  of  interaction  with 

the  bulk  GaAs  zone-center  LO  phonons,  both  in  magnitude  and  in  field 

dependence.  The  experiment  has  been  improved  and  extended  recently  in  a 

C  8  -T.\_ 

senes  of  measurements.  ihe  data  appear  to  deviate  sraoothlv  from  the 
calculated  transition  energy  in  the  absence  of  electron-phonon  interactions, 
and  the  deviation  is  generally  smaller  than  what  is  reported  in  Ref.  6.  In 


addition,  two  gaps  are  observed  and  manv  more  points  are  measured  above  the 

first  gap  in  qualitative  agreement  with  our  theoretical  results. 

we  report,  in  this  communication,  a  calculation  of  the  ls-2p  transition 

energy  as  a  function  of  the  applied  magnetic  field  for  a  magnetopoiaron  bound 

to  a  hydrogenic  impurity  in  a  double  heterostructre  CDHS)  i.e.,  a  single 

quantum  well  with  finite  barriers.  The  e 1 e c t r on - p honon  interaction 

Hamiltonian  derived  in  Refs.  4  and  5  is  treated  as  a  perturbation  on  the 

hydrogen- like  atom  confined  in  a  quantum  well  under  strong  magnetic  fields. 

Since  the  unperturbed  system  cannot  be  solved  exacity,  a  variational 

calculation  has  been  performed  with  properly  chosen  trial  wave  functions.  The 

calculation  is  actually  quite  complicated  and  we  can  only  present  the  results 

along  with  a  brief  outline  of  the  procedures .  Detailed  account  of  this  work 

9 

will  be  published  elsewhere. 

Consider  a  donor  impurity  at  the  center  of  the  GaAs  quantum  well  of 
width  d  in  a  GaAs/AiGaAs  DHS  system.  A  magnetic  field  is  applied  along  the 
growth  direction.  For  convenience,  we  introduce  the  two-dimensional  vectors  * 

“♦  “♦  "4  — ♦  — ¥ 

and  p  such  that  k  -  (*,q)  and  r  (p, z)  for  the  phonon  momentum  and  electron 

position,  respectively.  The  electron  momentum  is  denoted  bv  k  =  (k«,k  ).  In 

e  ,|  z 

the  absence  of  electron-phonon  interactions,  the  Hamiltonian  of  the  impurity 
is  given  by 


2 

e 

«  r 
o 


V3(z) 


VB(z) 


|  z I  <  d/2 
|  z |  >  d/2 


(2a) 


(2b) 
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The  interaction  Hamiltonian  H  ,  is  taken  directlv  from  Ref.  5  It 

e  -  p  h 

consists  of  two  terms:  the  electron  interaction  with  confined  LO  modes, 
and  the  interaction  with  interface  modes.  As  has  been  shown  in  Ref.  5, 
contributions  of  lattice  vibrations  outside  the  well  to  the  polaron  effect 
are  significant  only  when  the  well  width  d  is  extremely  small.  For  d  >  100  A, 
they  are  completely  negligible.  Therefore  we  have 


He-ph  "  He-L0  *  He-IN 


(3a) 


H  -  -  ■  e 
e  -  LO 


Lk'O  | 


3  (x)cos('"11*  z)[a  ik.)  -  a  -  x  j 
m  d  m '  m 


m-1 , 3 , . . . 


m—2 , 4 , 


B  (»c)sin(S7  z)[a  (*)  +  a‘(-x)]l 

m  dm  m  J 


(3b) 


H 


e-  IN 


\  aiK'P  Jr  /„■,  -F.9.sh..(k,zl_  , 

2  1  sj  cosh(xd/2)  1  s , j 


(")  +  “sj'-"'5 


„  /  x  sinh(xz)  ,  .-. 

B  .  (x)  .  ,  .  ,  '  a  .  (x)  +  a 

aj  sinh(xd/2)  a,j 


(3c) 


where  He  represents  the  electron  interaction  with  confined  bul.  LO  phonon 

modes  in  the  quantum  well  where  q  -  mw/d  is  quantized,  and  Hg  represents 

the  electron  interaction  with  interface  phonon  modes  in  the  well.  We  have 

introduced  in  Eqs .  (3)  the  creation  (annihilation)  operators  a(x)  (a  (x))  for 

m  ’  m 

>  *  -*  —* 
the  confined  modes  and  a1  .(x)  (a  .(x))  for  the  symmetric  and 

sj,aj  sj.aj 

antisymmetric  interface  phonon  modes.  They  obey  the  commutation  relations 


a  ( #c )  .  a  ,  (  k  '  '  '  =*  5  .  5  (k-k’  '■ 
'  at  o  20 


’  a '  (  k  )  ,  a  „  ( * '  '  '  -  a  <  * )  .  a  1  '  *  0 

a  j  a  3 


The  normalization  constants  are  given  by 


,  iire^trfuT  , 

B  (#c)  |  -  7“  — - “ - -  {-*-  -  — ' 

<  (mir/d)  *>1  01 


no)  .  i ,  k) 


B  .  .  k  )  |  ^  =  o  -  “  1  7” 

sj  .-vc  ;■ ,  cann  (  ko./  2 )  t-  e. 


„  2  )rfw  .  (*) 

d  fje,|2  .  a#-- - aj - - 

aj  Aye  £ 1  coth(/cd/2)  +  ' 


where  A  stands  for  the  interface  area,  e  and  e  denote  the  optic  and 

a>l/  oi/ 

static  dielectric  constants  of  material  v ,  0-1  for  the  well  and  i/  -  2  for 
the  barrier)  and  £  (w)  is  defined  bv 


£  (cj)  £  (w)  -  £ 

1/1/01/ 


:  ( w )  -  e 

i/  ®i/ 


with  the  dielectric  function  of  material  i/  given  by 


2  2  2  2 
«  (»)  -  e  Ot  )/(“t  ) 

i/  iii/  Li/  Ti/ 


we  now  proceed  to  calculate  the  transition  energy  bv  perturbation 


theory . 


For  the  unperturbed  system,  we  take  the  hvdrogenic  impurity  in  the 


6 


well  wich  a  magnetic  field  applied  normal  to  the  interfaces.  Thus  and 

H  is  treated  as  a  perturbation.  It  is  noted,  however,  that  even  the 

e  -  pn 

unperturbed  problem  is  not  exactly  solvable.  Thus  we  employ  the  variational 
method  to  determine  the  unperturbed  energy  levels.  The  trial  wave  function  is 
taken  to  be 


7 


£ . (B>  -  E. (B) 


2x)~ 

n»*i 


i 

d~  k 


£  .  (B)  -  E  (B)  -  ^(k) 
i  n 


(ID 


where  E.(B)  stands  for  the  corresponding  unperturbed  energy.  As  the  matrix 

element  is  generally  small  compared  to  the  transition  energy,  the  perturbation 

energy  in  (11)  becomes  appreciable  only  when  the  electronic  energy  level 

difference  in  the  denominator  matches  the  phonon  energy  hw(*) .  This  implies 

immediately  that  the  electron-phonon  coupling  has  negligible  influence  on  the 

Is  energy  level.  Hence  it  is  sufficient  to  calculate  AE  -  e „  -  E,  for  the 

2p+  Is 


transition  energy.  Furthermore,  we  note  that  significant  contribution  from 


the  second  term  in  (11)  is  expected  around  AE  =  E 


2p 


for  the 


three- level  resonance,  as  well  as  around  AE  =  for  the  two- level  resonance. 

As  can  be  seen  from  Eqs .  (3),  the  calculation  of  energy  levels  involves 
the  evaluation  of  four  matrix  elements 


•VT>  '  <2'V0JHe-Lolls'1»> 

•V7)  -  <2>V0jlH..INlls’1j> 

(12) 

*3<1)  -  <2p+.0jHe.L0|2p.,lm> 

H4(1)  -  <2p+,0.|He.I[1|2p..l.>  . 

where  the  states  Ini,  N  >  are  specified  by  the  atomic  level  ni  and  the  number 
'  m 

N  of  phonons  of  mode  m.  The  explicit  expressions  for  these  matrix  elements 


have  been  worked  out  in  Ref.  9.  Here  we  just  present  the  results  as  functions 
of  the  applied  field  in  Fig.  1.  It  is  observed  that  in  both  the  two-level  and 


3 


three-level  resonance  cases,  the  interface  phonon  modes  generally  result  in 
larger  matrix  elements,  ana  that  M.  is  the  largest  and  is  generally  2  to  a 

-4  ~ 

times  larger  than  the  others.  This  is  understandable ,  since  the  Bohr  radius 
of  the  impurity  is  of  the  same  order  as  the  well  width,  and  since  the  p-state 
wave  functions  extend  out  much  farther  in  the  xy-plane  than  the  s -state  wave 
function. 

The  energy  of  2p+  level  is  then  given  by 


«-,p  (7)  -  t.7)  +  K7) 


( 13a) 


T(  7) 


(  2jt)  ' 


d“« 


M, 


v€2p+  “Is  *"L  £, 


M, 


{2p  -Els*(E 


2p  •E1s^wL) 


£2p+'Els'(E2p/Els+^('c) 


d 


(13b) 


The  evaluation  of  (13)  is  still  very  difficult.  We  adopt  a 
iteration  procedure  to  compute  the  transition  energy.  Let  X( 7)  - 

and  T(7)  -  -  En ^ .  Then  Eq .  v!3)  can  be  rewritten  as 


nume  r i c  a  1 


X<7)  -  r(7)  +  T ( 7 , X ) 


(14) 


The  function  T(7)  represents  the  transition  energy  in  the  absence  of  the 
e lec  tron  -  phonon  interactions  and  is  already  known  from  the  variational 
calculations.  Therefore,  we  start  the  iteration  with  X  -  f(7)  in  T  for  a 
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fixed  B  which  determines  the  parameter  7.  A  new  X  value  is  oDtained  from  (Id) 

ana  employed  to  compute  a  new  I.  The  procedure  goes  on  until  seif  -  consistency 

is  achieved.  The  iteration  procedure  repeats  when  B  changes  and  eventually 

the  transition  energy  is  found  as  a  function  of  B. 

Results  of  our  numerical  computation  are  presented  in  Fig.  2  in  which 

the  is-2p^  transition  energy  is  plotted  as  a  function  of  B.  The  parameters 

used  in  the  numerical  work  are  d  -  125  A,  m  0.067  m  ,  a  -  98.7  A,  R  - 

e  00  y 

5.83  meV  and  V  -  230  meV.  The  solid  line  includes  the  effect  of  the 

o 

1 0 

nonparabolic  band  mass  and  the  dashed  dotted  line  does  not.  The  dashed  line 
represents  the  transition  energy  from  unperturbed  energy  levels.  It  is 
computed  from  Greene-Bajaj  wave  f  motions^  with  the  nonparaboiicity  effect 
included.  The  horizontal  lines  are  drawn  only  to  indicate  the  gap  positions. 
It  is  remarkable  that  the  theoretical  curve  breaks  into  three  branches 
separated  by  two  gaps,  in  good  qualitative  agreement  with  experimental  data. 
A  comparison  of  the  theory  with  experiments  can  be  found  in  Ref.  8.  The 
pinning  effect,  or  the  unperturbed  energy  level  repulsion,  is  apparently  a 
result  of  the  strong  resonance  interaction  when  the  denominator  of  the  second- 
order  perturbation  energy  vanishes.  As  the  electron  does  not  interact  with  TO 
phonons,  the  appearance  of  the  gaps  can  only  be  attributed  to  the  coupling  of 
the  electron  with  interface  phonons  which  oscillate  at  frequencies  between  the 
bulk  and 

Further  studies,  both  theoretical  and  experimental,  are  necessary  for  a 
complete  understanding  of  the  problem.  The  determination  of  the  transition 
frequency  from  experimental  data  is  difficult,  as  the  observed  intensity 
distribution  in  the  resonance  region  deviates  greatly  from  the  Lorentz  shape. 
Theoretical  study  of  the  influence  of  e 1 ec tr on  -  phonon  interactions  on  the 
transition  probability  is  needed  for  detailed  comparison  with  the  experimental 
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line  shapes  anc  observed  rapid  decrease  in  intensity  in  this  energy  region. 
Furthermore .  as  ve  have  noted  previously,  the  trial  function  11.  has  only  one 
adjustable  parameter.  More  accurate  functions  must  be  constructed,  for  more 
accurate  results. 

In  conclusion,  we  assert  that  it  is  the  interaction  of  the  electron  with 
interface  phonon  modes  that  change  the  polaronic  properties  in  the  reduced 
geometry.  These  modes  should  always  be  included  in  the  interpretation  of  such 
experiments.  Since  the  coupling  constant  of  the  bulk  Frohlicn  tvpe  cannot  be 
clearly  defined  for  the  interface  modes,  polaronic  phenomena  observed  in  the 
reduceo  dlmensionalitv  should  not  be  analyzed  bv  varying  the  couding  strength 
in  terms  of  the  bulk  Frohlicn  interaction. 

Ve  are  grateful  to  B.  D.  McCombe  and  J.  P.  Cheng  for  providing  us  with 
experimental  data  prior  tc  publication.  This  research  was  supported  in  part 
by  the  U.  S.  Office  of  Naval  Research. 
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Figure  captions 

Absolute  values  of  the  matrix  elements  of  the  electron-phonon  interactions 
calculated  as  functions  of  the  applied  magnetic  field  in  the  GaAs  well  of 
a  GaAs/GaAlAs  double  neterostructure . 

2.  ls-2p+  transition  energy  as  a  function  of  the  applied  field  B.  The  solid 
line  represents  results  from  the  present  theory  including  nonparablic 
band  mass,  and  the  dashed- do cted  line  does  not.  The  dashed  line  is  the 
transition  energy  in  the  absence  of  electron-phonon  interactions.  It  is 
calculated  from  Greene-Baiaj  wave  functions  with  the  nonparabolicity 
correction. 
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